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Abstract 

We study the oscillatory behaviour of a gene regulatory network 
with interlinked positive and negative feedback loop. Frequency and 
amplitude are two important properties of oscillation. Studied net¬ 
work produces two different modes of oscillation. In one mode (mode 
1) frequency remains constant over a wide range amplitude and in 
other mode (mode 2) the amplitude of oscillation remains constant 
over a wide range of frequency. Our study reproduces both features 
of oscillations in a single gene regulatory network and show that the 
negative plus positive feedback loops in gene regulatory network of¬ 
fer additional advantage. We identified the key parameters/variables 
responsible for different modes of oscillation. The network is flexible 
in switching between different modes by choosing appropriately the 
required parameters/variables. 
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1 Introduction 


The fundamental unit of life is the cell. Organisms may consist of just one 
cell or they may be multicellular. The multicellular organisms organized 
into tissues which are groups of similar cells arranged so as to perform a 
specihc function. One may view cell life as a collection of networks interact¬ 
ing through proteins, RNA, DNA and small molecules involved in signaling 
and energy transfer. These networks process environmental signals, induce 
cellular responses and execute internal events such as gene expression, thus 
allowing cells and entire organisms to perform their basic functions. These 
control and communication networks can be relatively simple (in bacteria) 
or they may be incredibly sophisticated (in higher organisms). In addition 
to their own needs for survival and reproduction, cells in multicellular organ¬ 
isms need additional levels of complexity in order to enable communication 
among cells and overall regulations. In living organism, proteins are the 
functional molecules. They are synthesized in a regulated processes known 
as Gene Expression (GE). So, gene expression and regulation are of fun¬ 
damental importance in cell. Again, proteins from one gene regulate the 
expression from other. In this way, gene regulatory networks have grown 
inside the cell. There can be other type of networks like metabolic networks, 
protein-protein interaction networks etc^. In general, the structure or ar¬ 
chitecture of the networks determines the function of the networks^. It is 
observed that positive and negative feedback loops are very common motifs 
in biological networks. They occur frequently in different gene regulatory and 
cell signaling circuits. In general, it is known that positive feedback loop in¬ 
duces a switch like behaviour and bistability and that negative feedback loop 
produces oscillations, suppresses noise/fluctuation effects etc. The loops are 
often coupled to perform various functions in the networks acting as bistable 
switches, oscillators, excitable devices etc.^“®. 

Rhythmic phenomena represent one of the most striking manifestations 
of dynamic behaviour in biological systems^°“^^. Gellular rhythms are gener¬ 
ated by complex interactions among genes, proteins and metabolites. They 
are used to control signaling, motility, growth, division and death. These 
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rhythms appear in many regulatory mechanisms that control the dynam¬ 
ics of living system. For example, neural and cardiac rhythms are associ¬ 
ated with the regulation of voltage dependent ion-channels, metabolic os¬ 
cillations originate from the regulation of enzyme activity and intracellular 
calcium oscillations involve the control of transport process while regulation 
of gene expression underlies circadian rhythms at the cellular level. There are 
some essential requirements for biochemical oscillations^^“^^. In the course of 
time, open systems that exchange matter and energy with their environment 
generally reach a stable steady state. However, once the system operates 
sufficiently far from equilibrium and when its kinetics acquires a sufficient 
nonlinearity, the steady state may become unstable. Feedback processes and 
cooperativity are two main sources of nonlinearity that favour the occurrence 
of instabilities in biological systems. When the steady state becomes unsta¬ 
ble, the system moves away from it, often bursting into sustained oscillations 
around the unstable steady state. Theoretical analysis shows that a nega¬ 
tive feedback networks with sufficient amount of time delay and nonlinearity 
produces oscillations^’^^’^^. The time delay in the networks can be created 
by a long chain of intermediate reactions or by an extra positive feedback 
loop. Different types of interlinked positive and negative feedback loops are 
observed in cellular systems with different number of nodes and links (Fig. 
1). Such coupled loops play a variety of roles, acting as bistable switches, os¬ 
cillators etc., although a single positive and negative feedback loop can also 
perform these functions under certain conditions. It is demonstrated that 
coupled or interlinked feedback loops are superior to single feedback loops 
as oscillators^®. A superior oscillator has the property of constant amplitude 
over a wide range of frequency. There may be another type of oscillation 
in which frequency remains constant though amplitude of oscillations may 
varied. Constant amplitude oscillations are important in heart beat, cell 
cycle etc. For circadian oscillations frequency should remain constant in dif¬ 
ferent environmental conditions. Here we study a gene regulatory network 
which show both type oscillations depending on the variation of appropriate 
parameter. 
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2 Interlinked gene transcription regulatory net¬ 
work: The model 


We consider a gene regulatory network consist of three genes X, Y and Z 
which synthesizes three proteins x, y and 2 : respectively. Three genes form 
a closed loop structure and the product of each gene represses the synthesis 
process from other in a cyclic way starting from X to Y to Z. In addition 
to that there is a autocatalytic positive feedback loop in X. The network 
architecture is identical to the module considered by Tsai et al.^®. Only 
difference is that, the network module considered by Tsai et ah is regulated 
at the degradation level but in our network the regulation is achieved at 
the synthesis level. The network is shown in Fig. 1. The dynamics of the 
network is driven by the following coupled nonlinear differential equation. 


dx j ^ ki I ^7 
dt~ ^ K2* + 


dt 


-k4y + 


ks 


( 2 ) 


dz j ^ k^ 

^ ^ + K^3 + yns 


( 3 ) 


The equations contain basically three different kinds of terms, viz., degra¬ 
dation, negative transcription or repression and autocatalysis. The oscilla¬ 
tory behaviour (Figs. 2 and 3) of the interlinked gene transcription regulatory 
network is studied by varying the different rate constants. We solve the cou¬ 
pled nonlinear equations to observe oscillation numerically by Runge-Kutta 
4 technique. To verify the stability of the network, we consider the random 
parameter values in the range given in Table 1. We observe that 500 out of 
8404 parameter sets (5.94%) yielded the oscillations in presence of positive 
feedback loop. But in absence of autocatalytic positive feedback loop we 
observe that 500 out of 7937 parameter sets (6.3%) yielded the oscillations. 
The last result (network without autocatalytic feedback loop) is completely 
different from the result of Tsai et al.^®. To study the role of positive feedback 
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loop in the network we measure frequency and amplitude of oscillations from 
the 500 oscillatory data sets with different rate constants as variable. Then 
we take a particular set of rate constant chosen from the 500 sets of param¬ 
eter values for which oscillations are observed. When we vary the repression 
strength (ki) on X, we observe that frequency remains constant over a wide 
range of amplitude. As the positive feedback strength increases the range of 
amplitude of oscillation over which the frequency remains constant increases 
(Fig. 4). But if we vary the degradation rate constant (k 2 ), we observe that 
amplitude remains constant over a wide range of frequency. The autocat- 
alytic loop in X increases that behaviour further (Fig. 5). Same behaviour 
is observed when both ki and k 2 varies simultaneously (Fig. 6). This obser¬ 
vation shows that the degradation rate has more impact on the oscillatory 
behaviour of the network. Fig. 7 shows that amplitude remains constant over 
a wide range of frequency when varied the autocatalytic positive feedback 
strength ky. 



Fig. 1. Gene transcription regulatory network with interlinked positive and 

negative feedback loop. 
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Fig. 2. Time dependent Oscillatory behaviour of X, Y and Z for the 
parameter values ki=266.152, k2=5.730, k3=331.660, k4=3.681, 
k5=494.232, k7=9.168, Ki=1.280, K2=0.982, K3=0.959, K4=18.882, 
ni=2.658, n2=2.048, n3=2.512, n4=3.940. 



Fig. 3. Different Phase plots corresponding to the oscillatory behaviour in 

Fig. 2 



Fig. 4. Amplitude versus Frequency plot when varied ki with different 
values of positive feedback strength shown in the graph. The other rate 
constants are hxed at A:2=5.730, /c3=331.660, A:4=3.681, A:5=494.232, 
A:7=9.168, 771=1.280, 772=0.982, 773=0.959, 774=18.882, ni=2.658, 
77-2=2.048, 77-3=2.512, n4=3.940. 
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Fig. 5. Frequency versus Amplitude plot with k 2 variation for different 
values of positive feedback strength shown in the graph. The other rate 
constants are hxed at A:i=266.152, A:3=331.660, /c4=3.681, A:5=494.232, 
A;7=9.168, Ari=1.280, ^2=0.982, ^3=0.959, ^4=18.882, ni=2.658, 
712=2.048, 713=2.512, n4=3.940. 



Fig. 6. Frequency versus Amplitude plot with /c 2 and /ci(both) variation for 
different values of positive feedback strength shown in the graph. The other 
rate constants are hxed at A:3=331.660, A;4=3.681, A:5=494.232, A;7=9.168, 
A:i=1.280, A:2=0.982, A:3=0.959, A:4=18.882, 771=2.658, 772=2.048, 

773=2.512, 774=3.940. 
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Fig. 7. Frequency versus Amplitude plot with A; 7 variation. The other rate 
constants are hxed at /ci=266.152, /i;2=5.730, /c3=331.660, /i:4=3.681, 
A;5=494.232, /c7=9.168, 771=1.280, 772=0.982, 773=0.959, 774=18.882, 
77-1=2.658, 772=2.048, 773=2.512, 774=3.940. 


Table 1; Different rate constants with their ranges used to solve the coupled 
equation numerically. 


Rate Constant 

Value/Range 

Description 

k2 

0-20 

Degradation rate constant 

k4^ 

0-20 

Degradation rate constant 

k& 

0-20 

Degradation rate constant 

ki 

0-500 

Negative feedback/repression Strength 

ks 

0-500 

Negative feedback/repression Strength 

h 

0-500 

Negative feedback/repression Strength 

^7 

0-100 

Autocatalytic Positive feedback strength 

77i 

2-4 

Hill coefficient of repression 

n2 

2-4 

Hill coefficient of repression 

ns 

2-4 

Hill coefficient of repression 

774 

2-4 

Hill coefficient of auto-activation 

Ki 

0-2 

Half maximum value of repressive Hill function at 2:=77i 

772 

0-2 

Half maximum value of repressive Hill function at x=772 

Ks 

0-2 

Half maximum value of repressive Hill function at y=Ks 

77i 

0-20 

Half maximum value of repressive Hill function at x=774 


3 Conclusion 

We study a gene regulatory network with interlinked positive and negative 
feedback loop. The loop produces two different modes of oscillation. In 























































one mode (mode 1) frequency remains constant over a wide range amplitude 
and in other mode (mode 2) the amplitude of oscillation remains constant 
over a wide range of frequency. For circadian rhythm, mode 1 oscillation is 
very important because organisms try to maintain a constant frequency of 
their daily clocks in spite of the variation of the environmental condition. 
Mode 2 oscillation is important for heart beat or cell cycle for which hxed 
amplitude of oscillations is very much crucial in different frequency region. 
Our study reproduces both features of oscillations in a single gene regula¬ 
tory network and show that the negative plus positive feedback loops in gene 
regulatory network offer additional advantage. We identihed the key param¬ 
eters/variables responsible for different modes of oscillation. The network is 
flexible in switching between different modes by choosing appropriately the 
required parameters/variables. Therefore, gene regulatory networks with in¬ 
terlinked positive and negative feedback loops work as more superior oscilla¬ 
tor rather than the signaling networks. 
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